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Abstract 

The heavy quarkonium spectrum of Two Color QCD (QC2D) at non-zero quark 
chemical potential fi and temperature T with fi/T ^ 1 has been calculated in both 
S- and P-wave channels using a lattice non-relativistic formulation of QC2D. As fi is 
varied, the quarkonium spectra reveal three separate regions, corroborating previous 
findings that there are three distinct physical regimes of QC2D at low temperature 
and high baryon density: hadronic matter, quark/quarkyonic matter, and deconfined 
matter. The results are interpreted in terms of the formation of heavy-light Qq states 
in the two-color baryonic medium. 
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1 Introduction 



Cold dense baryonic matter in Quantum ChromoDynamics (QCD) is difficult to study the- 
oretically because the lattice gauge theory method allowing a ffist principles investigation 
of non-perturbative physics does not work due to the "complex action problem": the intro- 
duction of a quark chemical potential /i 7^ in the Euclidean formulation of QCD makes its 
action complex, and importance sampling used for Monte Carlo evaluation of the partition 
function is difficult if not impossible to implement. 

In previous work [H |2l [3] we presented non-perturbative results at non-zero density 
for QC2D, a QCD-like theory with two colors, using orthodox lattice simulations. Since 
the gauge group is SU(2) and the quark representation pseudoreal, the functional measure 
det^^M(/i) remains real even once /i 7^ 0, and positivity can be ensured by insisting that 
the number of flavors Nf is even. In QC2D baryons are bosonic, and hadron multiplets, 
including Goldstone boson states associated with the breaking of global symmetries, contain 
both qq mesons and qq baryons. Despite these clear differences from physical QCD, QC2D 
has an unexpectedly rich structure as fi is increased, and we expect many lessons learned 
here about issues such as deconfinement, chiral symmetry restoration, and exotic ground 
states, may be more widely applicable. 

In brief, we have found three regimes with distinct behaviours beyond the onset 
at fi = fio = ""^77/2 where baryonic matter is ffist induced into the ground state at zero 
temperature. For fio ^ fJ' fJ'Q the matter consists of tightly bound qq scalars, which 
form a Bose-Einstein Condensate (BEC). Since the scalar diquarks are Goldstone modes 
associated with the spontaneous breaking of global baryon number conservation leading 
to superfluidity, thermodynamics in this regime is well-described by an effective approach 
based on chiral perturbation theory However, at the larger densities found in the range 
/^Q ^ ^ A^D, we find the thermodynamic properties of the system to be more like those of a 
degenerate system of quarks having a well-defined Fermi surface; superfluidity persists, but is 
now understood as resulting from a Bardeen-Cooper-Schrieffer (BCS) mechanism involving 
the condensation of weakly-bound diquark pairs at opposite points on the surface. Finally, 
for the largest densities /i > /id we find color deconfinement via the classic signal of the 
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non-vanishing expectation of the Polyakov loop. With the simulation parameters outlined in 
Sec. |2] below we estimate fio — 360MeV, /ig ~ 530MeV, and fio — 850MeV [2]. The regime 
f^Q ^ ^ fJ'D is particularly interesting, since it resembles "quark matter" while remaining 
color-confined; as such it is reminiscent of the "quarkyonic" phase first discussed in [5j. 

In this work, we study QC2D with both 7^ and temperature T > using corre- 
lation functions associated with heavy quarkonium QQ states, which offer a sophisticated 
probe of medium effects. The behaviour of quarkonia for T > in QCD is suggested as 
one of the signatures for quark-gluon plasma formation [B]. Since the energy required to 
create a. Q — Q pair is much larger than both Aqcd and the average thermal energy, the Q 
- Q production rate is dominated by short distance physics which is T-independent. Heavy 
quarkonium is a low energy QQ bound state, and the thermal medium found in, say, heavy- 
ion collisions can influence quarkonium formation. In [6], the authors modelled the thermal 
medium using a screening potential Vqq and showed the "melting" of heavy quarkonium 

(i.e., non-existence of QQ bound states above a certain temperature) by solving the result- 
ing Schrodinger equation. Evidence for suppression of excited bottomomium states at T > Tc 
in lattice simulations is given in [TJ 

It is important to stress that in this context the heavy quarks are regarded as test 
particles and are not in thermal equilibrium. In the current study this implies that Q, Q do 
not couple directly to the chemical potential fi. To treat the heavy quarks we will use an 
effective approach known as non-relativistic QCD (NRQCD); in QCD this is usually applied 
to the b - b system (the T and Xb meson families) where there is a well-deflned hierarchy of 
scales. For heavy quark mass M and velocity v, for NRQCD to be applicable we require [9] 

M >p ~ r"^ ~ > ~ Mw^ (1) 

where AE is the quarkonium "binding energy" and r is the typical distance between the 
quark and the antiquark. For T systems t;^ ~ 0.1. In Sec. [2] we briefly review the lattice 
approach to QC2D with /i 7^ and outline the formulation of NRQC2D. Our main results, 
for spin-singlet and spin-triplet states in both S- and P-waves are presented in Sec. El and 
a discussion of the observed /i- and T- variation follows in Sec. HI 
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2 Formulation 



We investigate the heavy quarkonium spectrum at non-zero temperature and baryon density 
by calculating 0{v'^) non-relativistic QC2D correlators using background lattice gauge field 
configurations generated on 16^ x 12, 12^ x 16, and 12^ x 24 lattices, at /3 = 1.9, k = 0.168 
with two dynamical fiavors of Wilson quark These parameters correspond to lattice 

spacing a = 0.186(8) fm (= 1/1.060(45) GeV"^), m^a = 0.68(1) and m^/nip = 0.80(1), 
where the scale is set by the string tension ((440 MeV)^ at /i = 0). The corresponding 
temperatures are T = 44, 66 and 88MeV. The range of chemical potential studied is < 
fia < 1.1; in |2] at /i = 0.8a^^ ~ 850MeV a quark density Uq = 16 - 32fm~^ was reported, 
corresponding to between 35 and 70 times matter density, where the uncertainty is due to 
discretisation artifacts. 

A standard Hybrid Monte-Carlo algorithm was used to generate lattice configurations, 
where the action 



M.,y = 5,,y - 5^ [(1 - 7.)e'^''"°f/.(a;)5,,.+. + (1 + ^,)e-^'-'>Uliy)Sy,,^,] . (3) 



The diquark source term proportional to j mitigates large infrared fiuctuations in a superfiuid 
phase with {^^2 i.^l^)'^'i''Pi) 7^ 0' ^^"^ ^1^° helps ergodicity by enabling real eigenvalues of M. to 
traverse the origin. To assess the effect of the diquark source, configurations generated with 
two different magnitudes j = 0.02 and 0.04 were used. Details of the simulation algorithm 
and previous analyses of various bulk thermodynamic quantities are given in [H [2l [3] and 
the temperature dependence of fin is discussed in [ini ITT] . 

We used the following 0{v^) Euclidean NRQC2D lagrangian density for the heavy 
quark with mass M: 



(2) 




with 



1/ 



(4) 
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with 



X, (5) 



and 



8M3 

+C2^ [0^(D-E-E-D)0 + xUD-E-E-D)x] 

-C3^ [0^ ■ (D X E - E X D) + ■ (D X E - E X D) x] 

-C4^ [0^^ ■ B0 - ■ Bx] (6) 

It is similar to NRQCD [12] with the only difference that and D are now gauge covariant 
temporal and spatial derivatives for SU(2) gauge theory. Here and x ^i-re two-color two- 
spinor fields for the heavy quark and anti-quark, and the tree- level value for the Cj is 1. 

We use the following discretised Green function of ([4]) to calculate the heavy quark 
Green function: 

G'(x,r = 0) = ^(x), 
G(x, r = a.) = - f/i(x, 0) (l " ^) " 0), 



G(x, r + = (l - t/i(x, r) (^1-^^(1- G(x, r), 

where 5'(x) is the source and the lowest-order hamiltonian reads (see fIlOp below) 



(7) 



and 



^0 = (8) 



5Hn = i- + . E - E ■ A) - -^o- • (A X E - E X A) 

8M3 8M2 ^ ^ 8M2 ^ ^ 



2M 24M 16nM2 ' 

5 



This approach is used for heavy quarkonium spectroscopy in QCD at T = in 
[T3| [HI [15] and recently in QCD at T > [7J. The integer n controls the high- momentum 
behaviour of the evolution equation. Since in QC2D we do not have any phenomenologically 
compelling choice for the heavy quark mass, we choose Ma > 3 and n = 1, to ensure 
compatibility with ([1]). The last two terms in SH are corrections to the kinetic energy term 
at non-zero lattice spacing [16]. The lattice covariant derivatives are defined by 

Ai0 = ^^Ui{x)(P{x + i)-Uj{x-i)(P{x-i) , 
A(2)0 = V = Y.^ \Ui{x)(P{x + i) - 2(P{x) + Uj{x- i)(P{x - i)] , 

i I 
i 

and E and B in Eq.Q are lattice cloverleaf definitions of the SU(2) chromoelectric and 
chromomagnetic fields. To mitigate quantum corrections we use tadpole improvement [T7] . 
replacing U^{x) Uq^U^{x), where Uq is the average link determined from the plaquette 
expectation value, and setting the coefficients q of ([6]) to 1. Note that in SU(2) gauge 
theory, heavy quarkonium (QQ) states are equivalent to heavy baryon states (QQ) [18] . 
Since the P-wave excitation needs an extended source. Coulomb gauge fixing is performed 
on SU(2) gauge fields prior to calculation of G(x, r). By combining non-relativistic heavy 
quark correlators, ^Sq, ^Si, ^Pq and ^Pi heavy quarkonium states could be studied. The 
expressions for the interpolating operators for all the states are listed in [19] . 

3 Results 

3.1 5- wave states 

We found that correlators for the S-wave states ^So and ^Si could be fitted with an expo- 
nential decay oc e~^^"^ corresponding to a simple pole even once /i 7^ 0; moreover the fits 
were quite stable over large ranges of r. This suggests that S-wave quarkonium bound states 
persist throughout the region < fj,a < 1.1 and ^ < Ta < j^. 
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Figure 1: (left) Energy of the ^5*0 state vs. quark chemical potential fi for heavy quark mass 
Ma = 3.0,4.0 and 5.0 with j = 0.02 on 12"^ x 24 lattice; (right) Temperature dependence of 
the ^5*0 state energy vs. fi for Ma = 5.0 with j = 0.04. 

In non-relativistic QCD, M„ = 2{ZmM — Eq) + AE'n for the state ra, where Zm is 
the heavy quark mass renormalization, Eq is a state-independent additive renormalisation, 
and AE'n is the fitted energy of the state [16]. Usually, the experimental value for one of the 
heavy quarkonium masses is chosen to fix Eq which is independent of n. This cannot be done 
in QC2D due to the lack of experimental spectrum data. However, since introducing 7^ 0, 
T > does not induce any new UV divergences, the change of the S'-wave state energy from 
that at /i = or T = 0, which must reflect underlying physics, can be measured. 

Fig. [T] shows the T- and //-dependences of the ^Sq state energy Ai?. The absolute 
values have unquantified contributions £'o(M), which in principle could be subtracted by 
matching AE{fi = 0), and as such contain little useful information. As argued above, 
however, the variation with fi is physical. As Ma increases, and hence higher order effects in 
v"^ become less important. Fig. [T] suggests three distinct regimes as fi is varied: initially the 
^5*0 state energy decreases from that at yU = 0, but once fi reaches the region /ii(— 0.5) < 
fia < yU2(— 0.85), the ^5*0 state energy stays roughly constant. For fi > fi2, the ""^^o state 
energy starts increasing again. The variation with /i becomes more marked as the heavy 
quark mass M is decreased; this may possibly be associated with the increasing size of the 
quarkonium state. The existence of three distinct /i-regions in which the S'-wave quarkonium 
state energies show markedly different behaviour and the agreement of fii^2 with the values 
fiQ and D found in [2] strengthens the argument for the existence of three different regimes 
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Figure 2: The splitting between the ^^i state energy and ^5*0 state energy for three different 
M on 12^ X 24 

described in Sec. [T]as the BEG phase, the BCS/quarlcyonic phase and the deconfined phase. 

The three temperatures shown in Fig. [1] all lie below the estimated deconfining tran- 
sition temperature T^a ~ | at /i = 0. As /i is increased, there is a deconfinement tran- 
sition, signalled by the Polyakov loop increasing from zero, at ^D{T)a ~ 0.3,0.55,0.75 
respectively [TOj [TT] . In accordance with this, we see that the ^5*0 en- 
ergy starts increasing roughly at the deconfinement transition at Ta = j^- However, 
no such behaviour is observed for the highest temperature. Interestingly, at fia = 0.6, 
AEa = 0.4157(6), 0.4226(2), 0.4235(2) as Ta rises from ^ to j^. This positive shift is simi- 
lar to that observed in the thermal mass of heavy quarkonium with increasing temperature 
in hot QCD with fx = 

The overall qualitative behavior of the ^5*1 state energy is quite similar to that of the 
^5*0 state energy. The ^Si state energy also shows the same three chemical potential regimes 
as the ^Sq state energy: the ^^i state energy decreases until fi reaches /ii and then stays 
roughly constant until /i reaches fi2 and then increases for fi > fi2- Thus, instead of the 



8 



0.46 




So a = 0.04) 
'S, a = 0.04) 



< 



0.41 







0.1 



0.2 



0.3 



0.4 



0.5 



0.6 



0.7 



0.8 



0.9 



1.1 



jLia 



Figure 3: diquark source dependence of the ^5*0 and the '^5'i state energies vs. fi with 



absolute ^5*1 state energy, the hyperfine sphtting AEas^ — AEis^ is shown as a function of /i 
in Fig. 121 Only a weak /i-dependence in the splitting for three explored Ms is observed, at 
most roughly 10% of the magnitude of the effect seen in Fig. [T] This may not be surprising 
since even in NRQCD at T = /i = 0, the hyperfine splitting is strongly affected by light 
quark dynamics and renormalization effects [21] , and the light quark mass in our simulation 
is relatively heavy {rriT^a = 0.68(1)). Clearly, much further study is needed to isolate chemical 
potential/temperature effects in the ^5*1 - ^5*0 splitting. 

The influence of the diquark source term on the heavy quarkonium spectrum and its 
effect in the three distinct /i-regimes we identified need to be studied. In general, increasing j 
induces a larger superfluid condensation pp. How this affects the heavy quarkonium spectrum 
is a highly non-trivial question. Thus, the spectrum calculation has been repeated with two 
different source magnitudes j = 0.02 and 0.04. Fig. [3] shows the j-dependence of AEisq and 
AEssj^- For both j = 0.02 and 0.04, the S'-wave state energies continue to manifest three 
separate regimes, but while for /i < /ii, the S'-wave state energies show little j-dependence, 
for fi > fii, the S-wave state energies for j = 0.04 are mostly larger than those for j = 0.02, 



Ma = 5.0 
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Figure 4: (left) The ratio Xlx ^(^' A*)/ Zlx '-^(^' 0) ^-^o correlators on 12^ x 24 with 
Ma = 5.0. Due to the noisiness of the P-wave data, only a limited r range is shown; (right) 
the corresponding ratio for ^5*0 correlators for comparison. 

signifying a larger diquark condensate effect in the quarkyonic and deconfined regions. 



3.2 P-wave states 

In contrast to the 5'-wave states, it is difficult to find stable exponential fits to the P-wave 
correlators with the current Monte-Carlo data before statistical noise sets in, except for 
the case fia < 0.25. We therefore use a more primitive quantity: the ratios of the ^Pq 
state correlators at several values /i 7^ (chosen within the three different regimes) to the 
correlator at /i = are compared in Fig. |H The other P-wave states (e.g., ^Pi state) show 
similar behaviors. The corresponding ratios of the ^Sq state correlators are shown to the 
right. Note that any effect we observe is entirely due to the dense medium. 

The S-wave correlators may be represented as a sum of exponential functions ^ ■ Aie~^ 
with AEi < AE2 < Ai?3 ■ ■ ■ so that the large-r behavior is dominated by the lowest energy 
state AEi. The S'-wave correlator ratios in Fig. H] confirm this expectation. At large r the 
ratio will be simply P(r; /i) = e'^'^f'^'' where (5(/i) = AEi(/i) - AEi{n = 0). In the BEC 
region (/i < /ii) and the deconfined region (/i > fi2), the ratio is ~ 20% at r/a = 23 and 
may be approximated as a straight line, R{t; /i) ~ 1 — 6{fi)T, consistent with the small, 
negative ^5*0 state energy difference 6{fi) that was previously observed. For the BCS region 
(/^i > /^fl > /^2)5 the ratio is ~ 60% at r/a = 23 and may be approximated as an exponential 
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function, which is consistent with the large energy difference in this region. In either case 
the ratio increases monotonically with r. 

Unhke the simple behavior seen for the S-wave, the P-wave correlator ratios show 
an interesting r- dependence. In the BEC and BCS regions {fi < ^2)1 the P-wave correlator 
ratios behave similar to the 5- wave, but in the deconfined region (yU > ^2)1 the P-wave 
correlator ratios are non-monotonic, initially decreasing with r before turning to rise above 
unity for ra ~ 4. On the other hand, the ratios of P-wave state correlators on 12^ x 16 and 
16^ X 12 lattices show monotonic behavior as the S'-wave correlator ratios do, which suggests 
a subtle interplay of density and temperature effects on the P-wave states. 

4 Discussion 

We have mapped out the variation of heavy quarkonium S'-wave states in a two-color baryon- 
rich medium with /i/T ^ 1, as both /i and T are varied. The behaviour is unexpectedly 
complex; the state energies initially decrease, then plateau, and finally rise again to become 
comparable or even exceed the vacuum value. The medium effect increases as the heavy 
quark mass M is decreased. Using the string tension to set the scale we find a downwards 
energy shift ~ 40MeV for Ma = 3.0 and yUi < /i < /i2. 

It is natural to seek an explanation in terms of the various ground states set out in 
[1112]. Initially assume T ^ 0; AE{fi) must then arise solely from interactions between the 
heavy quark pair QQ (or equivalently QQ) and light quarks q in the medium. An obvious 
possibility is the formation of two Qq states, or conceivably even a tetraquark QqqQ. In 
vacuum the quarkonium state usually lies (9(Aqcd) below the threshold for this to occur |22] . 
In a baryonic medium (ie. above onset /i > fio) the qs are already present and no longer 
need to be excited from the vacuum, so that the QQ energy may now be above threshold. 
A naive energy budget must take into account both the breaking apart of the gauge singlet 
QQ and qq states, and the subsequent formation of two Qq states. We assume that only qq 
breaking has any significant //-dependence, since it is the properties of the g-medium that 
evolve with /i. 

Henceforth assume that /^i coincides with the transition from BEC to quarkyonic 
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phase at hq identified in [2], and jj,2 with the deconfining transition aX hd- For fio < fi < fii 
the medium thus consists of tightly-bound diquark states, which are also Goldstone bosons 
associated with superfluidity, with mass proportional to ^yj [U [23]. We deduce that the 
energy required to break such bound pairs falls as the quark density Ug rises and, since 
the effect increases as M falls, with the ultimate separation of the resulting q - q system. 
Both these factors suggest an in-medium screening of the interaction between light quarks, 
or equivalently a non-trivial /i-dependence of the Goldstone decay constant F.„{fi), whose 
detailed mechanism remains unclear. 

For /ii < /i < the system is hypothesised to be in a quarkyonic phase, which 
we take to be a state in which quarks, though still confined, form a degenerate system 
with well-defined Fermi energy Ep ~ /i, and in which superfluidity arises through BCS 
condensation of weakly-bound and spatially delocalised Cooper pairs. To excite light qs 
capable of forming Qq states now thus requires an energy of 0{A), the superfiuid gap, 
believed to be approximately /i-independent in this regime [2]. Forming a Qq state at rest 
also requires the heavy quark to have kinetic energy 0{fi'^/M), which is a small, perhaps 
negligible, correction. This accounts for the approximate /i-independence of AE observed 
in this regime. Since the gs are no longer bound within Goldstone bosons, we expect their 
excitation energy in this regime to vary linearly with j, which is thus responsible for the 
mild increase of AE with j seen in Fig. [31 

Finally, for fi > fi2 the system is deconfined, as signalled by the non-vanishing expec- 
tation of the Polyakov loop. In this regime the physical states could in principle be isolated 
heavy quarks dressed by a cloud of both light quarks and now gluons, although we should 
not at this stage rule out the persistence of bound states - indeed, the contrast between S*- 
and P-wave states, which are more spatially extended, in this region in Fig. [H suggests an 
interesting story remains to be told. In either case the light constituents are now expected 
to have a Debye mass rriD o: gfi generated via quantum loop corrections; hence AE now 
rises with increasing /i, as confirmed by Fig. [TJ Thermal mass generation rriD oc gT may also 
be responsible for the systematic increase of AE with T seen in Fig. [T| it is notable in this 
case that the thermal effect appears to be equally manifest in all three /i-regimes. 

The reader will no doubt agree that these are speculative ideas, inevitably constrained 
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by theoretical pictures which can only become accurate in limits of vanishingly small or 
asymptotically high densities. It may well turn out that our identification of three different 
//-regimes is over-elaborate and masks a more unified explanation. Nonetheless, in view of 
the interesting heavy quarkonia physics in the non-zero temperature environment [3 [S], we 
expect from this exploratory study that heavy quarkonia can yield important insights into 
the nature of dense baryonic matter. 
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